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Abstract: Annular apodization of the illumination and/or imaging pupils of 
an adaptive optics scanning light ophthalmoscope (AOSLO) for improving 
transverse resolution was evaluated using three different normalized inner 
radii (0.26, 0.39 and 0.52). In vivo imaging of the human photoreceptor 
mosaic at 0.5 and 10° from fixation indicates that the use of an annular 
illumination pupil and a circular imaging pupil provides the most benefit of 
all configurations when using a one Airy disk diameter pinhole, in 
agreement with the paraxial confocal microscopy theory. Annular 
illumination pupils with 0.26 and 0.39 normalized inner radii performed 
best in terms of the narrowing of the autocorrelation central lobe (between 7 
and 12%), and the increase in manual and automated photoreceptor counts 
(8 to 20% more cones and 1 1 to 29% more rods). It was observed that the 
use of annular pupils with large inner radii can result in multi-modal cone 
photoreceptor intensity profiles. The effect of the annular masks on the 
average photoreceptor intensity is consistent with the Stiles-Crawford effect 
(SCE). This indicates that combinations of images of the same 
photoreceptors with different apodization configurations and/or annular 
masks can be used to distinguish cones from rods, even when the former 
have complex multi-modal intensity profiles. In addition to narrowing the 
point spread function transversally, the use of annular apodizing masks also 
elongates it axially, a fact that can be used for extending the depth of focus 
of techniques such as adaptive optics optical coherence tomography 
(AOOCT). Finally, the positive results from this work suggest that annular 
pupil apodization could be used in refractive or catadioptric adaptive optics 
ophthalmoscopes to mitigate undesired back-reflections. 
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1. Introduction 

Current clinical ophthalmoscopes form images of the retina by delivering and collecting light 
through non overlapping portions of the eye's pupil in order to avoid undesired reflections 
from the cornea and crystalline lens. Typically, the outermost part of the pupil is used for 
illumination and the central 2-3 mm are used for either observation (e.g. direct 
ophthalmoscope) or imaging (e.g. fundus cameras). When such a small area of the pupil is 
used for light collection, the image resolution is limited by diffraction when good spectacle 
correction has been achieved. For larger pupils, however, retinal images are blurred by the 
monochromatic aberrations of the optics of the eye. Ophthalmic instruments that incorporate 
adaptive optics (AO) [1-15] can compensate these aberrations, potentially reaching the 
classical Rayleigh resolution limit using larger diameter pupils. AO wavefront correction over 
a large pupil has resulted in an increase in transverse and axial resolution that has allowed in 
vivo visualization of multiple microscopic retinal structures in human and animal models, 
including but not limited to: the cone photoreceptor mosaic [1,6,16], the retinal pigment 
epithelial cell mosaic [8,17,18], leukocytes [19], the retinal capillary vasculature [19-21], 
fluorescently labeled retinal ganglion cells [8], the lamina cribrosa [22-24], and more 
recently, the rod photoreceptor mosaic [25,26]. Unfortunately, these structures cannot be 
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resolved in all subjects due to a number of factors, one of them being the reduction in pupil 
diameter that comes with aging [27] which results in a proportional reduction in transverse 
resolution and a quadratic reduction in signal. This issue alone severely limits the potential 
applicability of AO for clinical diagnosis, given that aging is a major risk factor in the most 
prevalent blinding conditions, including age-related macular degeneration, glaucoma and 
diabetic retinopathy [28-30]. It is therefore essential to explore methods for further improving 
on the transverse resolution provided by AO. 

Here, the use of pupil apodization, namely the manipulation of the amplitude and/or phase 
of the light at the pupil planes is explored as a means to narrow the point spread function 
(PSF) of the optical system. In particular, we evaluate the use of circular opaque masks 
centered on the pupils of the illumination and/or imaging paths [9,31,32] of a reflectance 
AOSLO, with the goal of improving its transverse resolution. The resolution improvement 
provided by these apodization masks measured in terms of the Rayleigh criterion [33], stems 
from the attenuation of the lower spatial frequency content of the images. Unlike other super- 
resolving techniques, such as structured illumination [34], pupil apodization does not increase 
the spatial frequency bandwidth of the image, and can in fact, be thought of as equivalent to 
edge-enhancing methods used in image processing. 

2. Theory 

The amplitude h\ of the PSF of a diffraction-limited rotationally symmetric lens with unit 
radius and a central obscuration of radius e under the paraxial approximation is given by [35] 



1 

h,(v,u)^^J,{vp)e- mp2l2 pdp, (1) 



where v and u are dimensionless radial and axial coordinates, defined as 

2k 



V = — 
X 



Here, A denotes light wavelength, / is the lens effective focal length, a is the lens radius (in 
this case also the mask outer radius), and x, y and z are the Cartesian coordinates with the z- 
axis along the lens axis of symmetry, which is also the direction of propagation of light. The 
combination of a rotationally symmetric lens with the type of mask described here is usually 
referred to as an annular lens. Figure 1 shows the normalized intensity of the PSF of such 
diffraction-limited annular lenses for the different values of e that will be evaluated 
experimentally. A small transverse narrowing and axial elongation of the PSF central lobe can 
be observed with increasing e. This means that any improvement in transverse resolution will 
come at the price of signal and axial resolution degradation. In confocal scanning instruments, 
the latter translates into a loss of axial sectioning, or equivalently, an increase in depth of 
focus. It should also be noted that this type of apodization increases the relative intensity of 
the PSF side lobes in the geometrical image plane. This should be kept in mind when imaging 
small dim objects in close proximity to brighter objects. 

A diffraction-limited confocal AOSLO is effectively a confocal scanning microscope, and 
as such, its PSF can be calculated by combining the effects of the imaging and illumination 
PSFs, as well as the finite size of the light detector. In an AOSLO, the effective size of the 
light detector is determined by the confocal aperture that precedes the light detector, while in 
AOOCT it is determined by the core of the optical fiber that is typically used to collect the 
light to be directed to the light detector. The intensity of the confocal microscope PSF is 
calculated as 



I OC |^j[i um i na tj on 
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Fig. 1 . Binary annular pupil transmission masks with various relative inner radii s (top) and the 
corresponding PSFs when used in diffraction-limited circularly symmetric lenses (bottom). The 
microscopy dimensionless coordinates u and v, defined in Eq. (2), are expressed in units of 
Airy disk radii (A.R). 

where p is a 2-dimensional function describing the detector size and <8> denotes convolution 
[36]. By combining Eqs. (1) and (3) we calculated the PSF for different values of e and 
detector size, as well as all four combinations of annular and circular lenses in the 
illumination and imaging pupils, and then determined the transverse and axial resolution in 
terms of the Rayleigh criterion [33]. The resulting resolution plots, shown in Fig. 2, do not 
take into account the wave-guiding properties of the photoreceptors that are used to explain 
the SCE, and therefore, departures from these theoretical predictions in our experimental 
results should be expected. 

2.1. Transverse resolution 

The most distinctive feature in all the plots in Fig. 2 is the transition that occurs when the 
detector size is comparable to that of the Airy disk. The resolution plots on the left remind us 
that in order to get the maximum improvement in transverse resolution provided by the use of 
a confocal aperture (up to 20%), its size has to be no larger than half the Airy disk. For AO 
ophthalmoscopes though, this benefit will only materialize if the AO compensation of the 
monochromatic aberrations is such that close to diffraction-limited performance is achieved. 
The transverse resolution benefit of using a sub-Airy disk confocal aperture comes at the price 
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Pupil masks Transverse resolution 



Axial resolution (a/f = 0.23) 



Detector signal 




Pinhole diameter (A.R.) Pinhole diameter (A.R.) Pinhole diameter (A.R.) 

Fig. 2. Transverse resolution (left), axial resolution (center) and energy detected (right) in a 
confocal scanning instrument, as a function of the detector size for various relative pupil 
obscuration radii (e), in terras of Airy disk radii (A.R.). The energy detected plots were 
generated assuming identical number of photons reaching the sample. 

of reduced light collection efficiency, as indicated by the plots on the right column, which 
could be a major limitation when the number of photons that can be delivered to the sample is 
limited. This is usually the case in fluorescence microscopy and in ophthalmoscopy where 
bleaching and thermal or photochemical damage to the retina, are a concern. 

The comparison of the transverse resolution curves for apodizing masks with different e 
shows that the apodization of the illumination path improves resolution for all detector sizes, 
with the largest benefit taking place for the larger confocal aperture. Interestingly, apodizing 
the pupil in the imaging path provides no benefit other than for confocal apertures much 
smaller than the Airy disk. The asymmetry between apodizing the illumination or the imaging 
pupils stems from the convolution of the imaging PSF and the detector function, shown in Eq. 
(3). The convolution effectively blurs any narrowing of the PSF of the apodized imaging 
optics other than for detectors smaller than a quarter of an Airy disk, for which the resolution 
improvement is approximately 5%. But in obtaining this small benefit the detected signal is 
reduced by almost an order of magnitude, which could lead to an unacceptably low signal-to- 
noise ratio (SNR). The plots that correspond to the simultaneous apodization of both the 
illumination and imaging pupils are similar to those that result from apodizing the 
illumination pupil only, with the added disadvantage of the detected signal reduction that 
comes with the annular imaging pupil. The sole situation in which the apodization of both 
pupils shows an improvement over the apodization of the illumination pupil alone (3% 
for e — 0.52 ), is for confocal apertures smaller than one quarter of the Airy disk. 
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2.2. Axial resolution 



The main advantage of using confocal scanning instruments for retinal imaging comes from 
the axial sectioning that fundus cameras cannot provide. The finite size of the confocal 
aperture in the image plane means that some out-of-focus light is not detected, thus enhancing 
the contrast of the in-focus feature. The plots in the middle column of Fig. 2 show that axial 
resolution for all values of e and pupil apodization configurations is better for detectors 
comparable to or smaller than the Airy disk. Unlike the transverse case, there is no axial 
resolution benefit in using confocal apertures much smaller than the Airy disk. 

Consistent with the PSF axial widening due to increasing e in annular lenses (shown in 
Fig. 1), the apodization of the illumination pupil in a confocal scanning instrument results in 
degradation the of axial resolution. This increase of the depth of focus could be beneficial for 
applications such as AOOCT. Apodization of the imaging path, however, does not affect the 
axial sectioning, other than for detector sizes smaller than the Airy disk. Finally, when both 
pupils are apodized, the degradation in axial resolution is approximately the sum total of the 
resolution loss that results from apodizing each pupil separately. 

2.3. Modulation transfer function 

The origin of the changes in transverse resolution with the inner radius and configuration of 
the pupil apodization mask can be interpreted in terms of its effect on the spatial frequency 
content of the image, as illustrated by the plots in Fig. 3. Typically, MTF curves are 



Pupil masks 
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□ 

Illumination 
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Imaging 

□ 
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Fig. 3. MTF of a confocal microscope with different apodizing binary mask combinations (one 
per row) with different radii (indicated by color) and pinhole sizes (one per column). These 
plots were calculated assuming 680 nm light, an eye with 17 mm focal length imaged through a 
7.75 mm pupil diameter. The cut off spatial frequency (f mm ) for an incoherent imaging system 
with uniform illumination (i.e. a flood illumination fundus camera) is provided for reference. 
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normalized to unit amplitude at zero frequency. In the plots presented here, however, the 
amplitude at zero frequency is proportional to the detected power, assuming that the same 
number of photons per second is delivered to the sample. In this way, we facilitate the 
comparison of MTF amplitude for different apodization configurations at any given 
frequency, given the same light exposure. The curves were generated assuming that 680 nm 
light was used to image an eye with a 17 mm focal length through a 7.75 mm pupil diameter. 

The plots show that for confocal pinholes smaller than the Airy disk the MTF extends 
beyond the cut off spatial frequency (/max) of a non-scanning incoherent imaging system that 
uses uniform illumination at the object plane (e.g. AO fundus camera). This increased 
bandwidth is the reason for the superior lateral resolution of confocal microscopes over bright 
field microscopes. When using detectors equal to or larger than the Airy disk both 
microscopes have the same theoretical transverse resolution, although in practice the confocal 
microscope is preferred because of its axial sectioning. The MTF curves in Fig. 3 also show 
that the use of an annular illumination pupil enhances the higher spatial frequency content 
slightly while maintaining or attenuating the intermediate and lower spatial frequencies. It is 
worth noting that the range of frequencies that are attenuated is between 0 and /max (1 + s)/2, 
thus indicating that the value of e used should be chosen with the spatial frequency content of 
the object of interest in mind. 

The use of annular imaging pupils reduces the MTF amplitude across its whole bandwidth 
but more so for the higher spatial frequency range, when compared to the MTF of the 
unapodized system. This is consistent with the fact that no transverse resolution improvement 
is predicted for this apodization configuration (see second row, leftmost plot in Fig. 2). 

The changes in the MTF of the simultaneously apodized illumination and imaging pupils 
can be coarsely described as those introduced by the annular illumination attenuated by the 
annular imaging pupil. 

3. Methods 

Images of the photoreceptor mosaic were collected in two subjects with no known history of 
eye disease. Written informed consent was obtained after the nature and possible risks of the 
experiment were explained to the subjects. The study was approved by the Institutional 
Review Boards at the University of Rochester and the Medical College of Wisconsin. The 
pupil of the imaged eye was dilated and cycloplegia was induced by using one drop of a 
combination of phenylephrine hydrochloride (2.5%) and tropicamide (1%), applied topically. 
The subjects were aligned and stabilized in front of the AOSLO using a dental impression on 
a bite bar. Photoreceptor images were collected at the fovea and at a rod dominated region 
approximately 10° temporal to the foveal center (determined by fixation). 

The two identical AOSLOs used in this study were described previously [10]. The only 
modification was the placement of circular opaque masks centered in pupil conjugate planes 
of the illumination (immediate to the light source) and/or imaging paths (immediate to 
detector). The illumination and imaging paths of the wavefront sensor were not apodized, 
using full circular pupils both in the illumination and imaging paths. The mask's central 
obscuration was held by a single, 0.5 mm thick post. Masks with relative radii e = 0.26, 0.39 
and 0.52 were used. Whenever the illumination and imaging pupils were apodized 
simultaneously, equal sized masks were used. Images for the different apodization 
configurations were collected within minutes of each other for a given e, in random order. 

The light collected through a 1.0 Airy disk diameter confocal aperture was detected by a 
photomultiplier tube (PMT) with gain adjusted to maintain an approximately constant mean 
pixel value across apodization configurations. An additional set of image sequences were 
collected at a fixed PMT gain to analyze the reduction of photoreceptor reflectance signal due 
to the SCE. A sup erlumine scent diode with a peak wavelength of 680 nm and 8.5 nm full- 
width at half-maximum bandwidth (46 urn estimated coherence length) was used as the 
imaging light source. This source delivered a fixed power to the illumination pupil that in the 



#164634 - $15.00 USD Received 13 Mar 2012; revised 25 Apr 2012; accepted 1 Jun 2012; published 20 Jun 2012 
(C) 2012 OSA 1 July 2012 / Vol. 3, No. 7 / BIOMEDICAL OPTICS EXPRESS 1654 



absence of an obscuration mask was 150 |iW at the eye over a 7.75 mm diameter pupil. The 
images were recorded over a 1.0° square field of view. At every retinal location and for each 
s, sequences of 200 images were recorded at approximately 17 frames per second for each of 
the four possible pupil configurations. This amounted to a maximum retinal radiant exposure 
five times below the ANSI maximum permissible exposure [37]. Fifty images from each 
sequence were registered using the method outlined in [38] and then averaged to increase 
SNR. 

4. Results and discussion 

4.1. Foveal photoreceptor mosaic 

Registered averages of image sequences captured at 0.5° from the foveal center are shown in 
Fig. 4 where some minor resolution improvement due to the use of apodizing masks in the 
illumination path can be observed for e = 0.26 and 0.39. Careful examination of the image 
corresponding to £ = 0.52 and annular pupil illumination (rightmost image on the second row) 
reveals that the intensity profiles of the cone photoreceptors appear multi-modal. These multi- 
modal patterns were first observed ex vivo by Enoch [39] and in vivo in humans by Dubra et 
al. [25] and are thought to be due to the optical wave-guiding properties of photoreceptors. 

Attempts to quantify changes in transverse resolution or image contrast using global 
sharpness metrics, such as those proposed by Muller and Buffington [40], proved inconclusive 
due to the temporal fluctuations of the photoreceptor reflectance [11,41,42]. A coarse 
estimation of the size of photoreceptors can be obtained by measuring the width of the central 
lobe of the image normalized autocorrelation, using the normalized correlation of the image 
with itself as defined in Ref. [38]. The change in full-width half-maximum (FWHM) values of 
the central lobe of the autocorrelation as a result of using annular pupils for all the images in 
Fig. 4 are shown in Table 1 below. These numbers loosely agree with the theoretical 
predictions that apodization of the illumination pupil translates to a sharpening of the cone 
photoreceptor intensity profiles and apodization of the imaging pupil does not. The low signal 
in the raw images that result from apodizing both pupils degrade the quality of the 
registration, resulting in blurrier images, as is most noticeable for e = 0.39 and 0.52. 

As predicted by the MTF calculations, the image radial spectra and their ratios (shown in 
Fig. 5) show an enhancement of the higher spatial frequencies below the cutoff frequency f mdx 
relative to the lower spatial frequencies. This enhancement is negligible for the apodization of 
the imaging pupil and maximal for the apodization of the illumination pupil, as predicted by 
the theoretical calculations. Interestingly, because the range of (low) frequencies that are 
attenuated increases with e, choosing a mask with too large a central obscuration can have the 
detrimental effect of attenuating the spatial frequencies that correspond to the object of 
interest (in this case, the photoreceptors), as it appears to be the case for e = 0.52. 

As expected, obscuring the center of either pupil affects the amount of light reaching the 
detector. To quantify this change and compare it to the values predicted by the SCE, the 
imaging was repeated with constant PMT gain. Despite the photoreceptor reflectance 
fluctuations [11,41,42] and other sources of variability such as image registration, tear film 
and subject alignment, the observed reflectance amplitude values (shown in Table 2) are 
consistent with the SCE predictions. 

The practical use of the proposed apodization for in vivo imaging of the photoreceptor 
mosaic should ultimately be assessed based on whether it improves our ability to identify 
photoreceptors. To this effect, we performed cell counts using a semi-automated algorithm 
developed by Li and Roorda [44] and further automated by Garrioch et al. [45]. Manual 
counting was also performed for comparison, twice on every image by three trained observers. 
The results from both counting methods, summarized in Table 3, are in agreement with the 
previous observations. The use of apodization masks with e = 0.26 and 0.39 in the 
illumination path allows the identification of 8 to 20% more cone photoreceptors, with both 
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Fig. 4. Registered averages of foveal cone photoreceptor AOSLO images. Each column 
corresponds to a different apodization mask normalized radius £, while each row corresponds to 
a different apodization configuration. For each £ all images show the same retinal location. 
Image contrast has been stretched for display purposes only. Scale bar is lO^m across. 

Table 1. Normalized autocorrelation central lobe FWHM changes due to pupil 
apodization corresponding to images in Fig. 4. Negative values indicate FWHM reduction 
that can be interpreted as sharpening of the image features. 



Pupil masks 


£ = 0.26 (%) 


£ = 0.39 (%) 


£ = 0.52 (%) 


None 


0.00 


0.00 


0.00 


Illumination 


-7.1 


-12.0 


-4.9 


Imaging 


3.0 


1.0 


4.6 


Both 


-7.1 


-1.2 


20.0 
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S = 0.26 



8 = 0.39 



S = 0.52 




Fig. 5 . The top row shows the radial spectra of the foveal cone photoreceptor images in Fig. 4, 
with each column corresponding to a different apodizing mask inner radius s, normalized to 
unit amplitude at zero frequency. The second row shows the ratios of the normalized spectra 
from images acquired with annular pupil(s) to the ones with full circular pupils. The cutoff 
spatial frequency for a non-scanning incoherent imaging system and the peak frequency 
associated to the cones are indicated for reference. 



Table 2. Measured relative average image intensities for images collected with constant 
PMT gain compared against predicted values (in parenthesis) according to the SCE 

measured objectively at fixation (p = 0.11 mm" 2 ) [43]. A Gaussian intensity profile at the 
illumination pupil (o = 2.9 mm) and a one Airy disk confocal aperture were assumed. 



Pupil masks 


£ = 0.26 


£ = 0.39 


£ = 0.52 


None 


1.00(1.00) 


1.00(1.00) 


1.00(1.00) 


Illumination 


0.72 (0.73) 


0.50 (0.49) 


0.21 (0.27) 


Imaging 


0.43 (0.64) 


0.27 (0.37) 


0.27(0.17) 


Both 


0.38 (0.47) 


0.17(0.18) 


0.03 (0.05) 



Table 3. Cone counts for the images in Fig. 4. The manual counts are averages (and 
standard deviation) from 6 counts by 3 different observers. 



Pupil masks 


£ 


= 0.26 


£ 


= 0.39 


£ 


= 0.52 


Automated Manual 


Automated Manual 


Automated Manual 


None 


283 


286 (1.3) 


266 


270(1.5) 


253 


266(2.1) 


Illumination 


306 


309(1.6) 


318 


322 (2.2) 


380 


385 (2.7) 


Imaging 


283 


287 (2.3) 


268 


273 (1.3) 


270 


276 (2.6) 


Both 


301 


306 (1.5) 


323 


329(1.8) 


350 


354(1.9) 



the automated and manual counting. The multi-modal appearance of the cones for the image 
corresponding to an annular illumination pupil with e = 0.52 led to a large number of false 
counts, highlighting the fact that the inner radius of the apodizing mask must be selected with 
some care. The apodization of the imaging pupil with the two smaller masks yields 
comparable cone counts to those resulting from the unapodized configuration, as expected 
from the theoretical predictions. 
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Interestingly, the multi-modal intensity patterns of the cone photoreceptors are either not 
as dramatic or not observed at all when using a circular illumination pupil and an annular 
imaging pupil. This is probably not due to an asymmetry in the SCE, but to the blur induced 
by the convolution of the imaging PSF and the confocal pinhole transmission in Eq. (3). 

4.2. Photoreceptor mosaic at 10° temporal 

At 10° temporal from fixation (Fig. 6) the rod photoreceptors are much smaller than the cones. 
As with the foveal images, more photoreceptors can be identified when the illumination pupil 

Pupil masks S = 0.26 S = 0.39 8 = 0.52 



Illumination 

□ 

Imaging 

□ 




Illumination 

□ 

Imaging 

□ 

Illumination 

□ 

Imaging 

□ 

Illumination 

□ 

Imaging 

□ 




Fig. 6. Registered averages of photoreceptor AOSLO images at 10° temporal from fixation. 
Each column corresponds to a different apodization mask normalized radius £, while each row 
corresponds to a different apodization configuration. For each £ all images show the same 
retinal location. Image contrast has been stretched for display purposes only. Scale bar is lOum 
across. 



#164634 - $15.00 USD Received 13 Mar 2012; revised 25 Apr 2012; accepted 1 Jun 2012; published 20 Jun 2012 
(C) 2012 OSA 1 July 2012 / Vol. 3, No. 7 / BIOMEDICAL OPTICS EXPRESS 1658 



is apodized. The radially averaged power spectra of these retinal images and their 
ratios,shown in Fig. 7, closely resemble those in Fig. 5 and are consistent with the theoretical 
calculations. 

The changes in autocorrelation central lobe FWHM as a result of the pupil apodization 
(shown in Table 4) are consistent with those observed when imaging the foveal cones. 
Interestingly, the largest change occurs when using the e = 0.39 mask which was also the case 
for foveal cones, suggesting that this might be the best of the three masks (and configuration) 
for photoreceptor imaging. 

Previous SCE studies show both through objective and subjective measurements that 
peripheral cone photoreceptors are more directionally sensitive than cones near the center of 
the fovea [39,46], and that the directional sensitivity of rod photoreceptors is similar to that of 
foveal cones [47,48]. By repeating the imaging shown in Fig. 6 with constant PMT gain, we 
recorded the mean relative photoreceptor intensity values presented in Table 5 below, that 
suggest an SCE described by lower p values than those published in the literature. This could 
be explained by the fact that the AOSLO images are focused on the reflections from the inner- 
outer segment interface, which according to the work by Gao et ah [49] are less directional 
than those from the posterior surface of the photoreceptor outer segments. As a consequence 
of the difference in SCE between rods and cones, the ratio of cone to rod signal decreases as e 
increases (see Table 6). This relative intensity change could be exploited to automate the 
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Fig. 7. Radial spectra of the rod-dominated images in Fig. 6, with each column representing a 
different £. The second row shows ratios of the normalized spectra from images taken with 
annular pupil(s) to the spectra of images taken with circular pupils. 



Table 4. Normalized autocorrelation central lobe FWHM changes due to pupil 
apodization corresponding to images in Fig. 6. Negative values indicate FWHM 
reduction, that can be interpreted as sharpening of the image features. 



Pupil masks 


£ = 0.26 (%) 


£ = 0.39 (%) 


£ = 0.52 (%) 


None 


0.00 


0.00 


0.00 


Illumination 


-8.0 


-12.7 


-7.4 


Imaging 


-3.7 


-2.2 


3.3 


Both 


-11.2 


-10.5 


3.1 
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Table 5. Measured relative average photoreceptor intensities for images collected with 
constant PMT gain compared against predicted values (in parenthesis) corresponding to 
Fig. 6. The values used to describe the SCE are the best fit to our experimental data (p cme! 
= 0.15 mm" 2 and />rods — 0.11 in n 1 ) assuming a Gaussian intensity profile at the 
illumination pupil (<r = 2.9 mm) and a one Airy disk confocal aperture. 





e = 


0.26 


e = 


0.39 


e = 


0.52 


Pupil masks 


Cones 


Rods 


Cones 


Rods 


Cones 


Rods 


None 

Illumination 

Imaging 

Both 


1.00(1.00) 
0.48 (0.68) 
0.59 (0.59) 
0.26 (0.40) 


1.00(1.00) 
0.55 (0.73) 
0.72 (0.64) 
0.35 (0.47) 


1.00(1.00) 
0.46 (0.42) 
0.33 (0.30) 
0.15(0.13) 


1.00(1.00) 
0.57 (0.49) 
0.36 (0.37) 
0.22(0.18) 


1.00(1.00) 
0.15(0.21) 
0.11 (0.12) 
0.02 (0.03) 


1.00(1.00) 
0.20 (0.27) 
0.13 (0.17) 
0.04 (0.05) 



Table 6. Ratios of the average image intensity from cones to that of rods 



Pupil masks 


s = 0.26 


£ = 0.39 


£ = 0.52 


None 


1.00 


1.00 


1.00 


Illumination 


0.87 


0.81 


0.75 


Imaging 


0.82 


0.92 


0.84 


Both 


0.74 


0.68 


0.50 



identification of cones and rods as follows. Two images one with and one without apodization 
can be collected and once each cell has been counted, the decision of which cell is a cone or a 
rod can be made based on the relative change of brightness between the two images. 

Lastly, manual cone counts, followed by automated and manual rod counts were 
performed on the images in Fig. 6. Once the cones were manually counted, the automated 
software counted all the bright spots as individual photoreceptors. Then the counts within 6 
pixels of the manually entered cone locations were removed as misidentified rods. The 
remaining identified bright spots were counted as rods. The results, summarized in Table 7 
further support the previous observations, showing the highest rod count for the imaging 
configurations with apodized illumination only, followed by the case with both illumination 
and imaging apodized. No consistent benefit or degradation was observed from apodizing 
only the imaging pupil. As with the cone images, the increase in photoreceptor counts is non- 
negligible (between 10 and 29%) when using annular illumination pupils, demonstrating the 
practical utility of the method. 

Table 7. Rod counts for the images in Fig. 6. The manual counts are averages (and 
standard deviation) from 4 values by 2 different individuals. 

Pu il masks c = 0-26, 41 cones £ = 0.39, 43 cones c = 0.52, 50 cones 

Automated Manual Automated Manual Automated Manual 

None 428 441 (4) 412 426 (5) 415 428 (6) 

Illumination 477 489(4) 533 549(10) 499 507(2) 

Imaging 393 404(5) 439 458 (8) 397 408 (4) 

Both 452 461 (3) 523 530 (2) 508 515 (6) 

5. Conclusions 

The use of centered, annular pupil illumination and/or imaging pupils was evaluated in a 
reflectance AOSLO as a method for increasing transverse resolution with low cost and 
complexity. Images from the photoreceptor mosaic at retinal eccentricities with very different 
spatial frequency content were collected and analyzed. Despite not considering the wave- 
guiding or scattering properties of the photoreceptors in the theoretical calculation, the 
experimental findings are in agreement with the predictions from paraxial confocal 
microscopy theory. This indicates that the use of annular pupil illumination and full circular 
imaging pupil provide highest transverse resolution. The two annular illumination pupils with 
smaller inner radii performed best, whether measured in terms of the narrowing of the 
autocorrelation central lobe or the manual and automated photoreceptor counts. The use of the 
annular apodizing mask with the largest inner radius (e = 0.52) provided different outcomes 
depending on the size of the photoreceptors; when imaging rods, the use of the mask 
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improved their visibility as measured with the autocorrelation and counts, while when 
imaging cones their intensity profile could become multi-modal. These multi-modal intensity 
profiles could lead to incorrect high cell counts as each lobe within the pattern might be 
identified as a rod or another cone. Therefore, the size of the apodizing mask must be chosen 
with some knowledge of the spatial frequency content and optical properties of the structure 
to be imaged. 

The appearance or enhancement of multi-modal patterns when illuminating the pupil with 
a non-uniform intensity profile is consistent with the wave-guiding properties of 
photoreceptors and it could be exploited to gain insight into the structure of photoreceptors 
and their dynamics. For example, by placing a spatial light modulator in the pupil plane and 
collecting images with different intensity patterns (not just annular, as demonstrated in this 
work), information on the modes of each photoreceptor could be obtained and used in 
combination with wave-guiding theory to infer the geometry of each photoreceptor. 
Alternatively, a fixed illumination pattern (e.g. annular) could be used to monitor the 
evolution of the multi-modal patterns, to determine if subtle transverse changes take place in 
the photoreceptor mosaic. This could be used to study the disc shedding and renewal process, 
in a similar way that the evolution of the fundamental modes has been used before [1 1,42,50— 
52]. If the multi -modal intensity profiles are not affected by daily changes, they could 
potentially reveal early changes in photoreceptor disease such as swelling. 

We found that the use of apodizing masks affected the average photoreceptor intensity as 
predicted by the SCE measurements reported in the literature. The fact that the SCE is 
different for cones and rods at any given retinal location means that cone intensity can be 
differentially attenuated with respect to that of rods through pupil apodization. By capturing 
pairs of images of the same retinal location with and without pupil apodization, the cone 
photoreceptors could be differentiated from rods. One potential imaging protocol is as 
follows: first, capture an image sequence with annular pupil illumination and e between 0.3 
and 0.4 for best transverse resolution. Next capture a pair of image sequences with and 
without apodization in both illumination and imaging pupils. A difference map, created from 
the pair of sequences can be used to identify the photoreceptors in the first image sequence. 

A consequence of the use of the apodization masks evaluated in this work is an axial 
elongation of the PSF. Although this fact can be seen as detrimental in confocal scanning 
ophthalmoscopy, it could be beneficial for AOOCT. The axial resolution of most current OCT 
ophthalmoscopes is proportional to the bandwidth of the source and mostly independent of the 
axial PSF, with current instruments boasting axial resolution of up to 3 \im. When imaging the 
retina using AOOCT with a full pupil, the lateral resolution degrades axially outside the PSF 
[13]. The use of annular illumination would extend the PSF axially and thus increase the 
range of retinal layers that can be imaged simultaneously without loss of transverse and axial 
resolution. In fact, the former would be improved. 

Finally, undesired back-reflections from the surfaces of refractive elements in most current 
AO ophthalmoscopes are removed by slightly tilting the elements. Placing an obscuration in 
the illumination path as demonstrated here could be used to mitigate these back-reflections 
without having to tilt these elements [53], thus potentially improving the optical performance 
of refractive AO ophthalmoscopes. The on-axis nature of refractive systems, coupled with the 
additional surfaces for aberration control might allow for smaller instruments. 
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